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, Abstract. For the first time an excess of photons above an energy threshold of 730 GeV from the giant radio galaxy M 87 

i-S^ ' has been measured at a significance level above 4cr. The data have been taken during the years 1998 and 1999 with the 

HEGRA stereoscopic system of 5 imaging atmospheric Cherenkov telescopes. The excess of 107.4±26.8 events above 730 GeV 
Q . corresponds to an integral flux of 3.3 % of the Crab flux or Ny(E > 730 GeV) = (0.96 ± 0.23) x 10"'^ phot. cm"-s"'. M 87 is 

( ' located at the center of the Virgo cluster of galaxies at a relatively small redshift of z = 0.00436 and is a promising candidate 

, among the class of giant radio galaxies for the emission of TeV y-radiation. The detection of TeV y-rays from M 87 - if 

■ confirmed - would establish a new class of extragalactic source in this energy regime since all other AGN detected to date at 

TeV energies are BL Lac type objects. 
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1. Introduction lAharonian et al. 2002a> have also been established as TeV y- 

ray emitters. However, other types of AGN, e. g. giant ra- 
Active Galactic Nuclei (AGN) are believed to contain as galaxies, also show relativistic mass outflows, though, in 
a central "engine" a supermassive black hole which causes ^^^^^^^^ Lac type objects, under large viewing angles, 
the development of large scale jets. Extragalactic TeV y-ray Amongst fliese the nearby radio galaxy M 87 has been spec- 
emission has been observed so far from AGN only of the ^^j^^^^ ^ powerful accelerator of cosmic rays (includ- 
BL Lac type, i.e. objects ejecting matter in a jet oriented j„g highest energy particles observed in the universe, see 
very close to the observer's line of sight. In BL Lacs, TeV ^ g Qinzburg & Syrovat skii (1964) Biermann et al. (2000)}. 
photons are commonly believed to originate in the relativis- ^ 87 has been targeted with the HEGRA Cherenkov telescopes 
tic jets, most popularly due to inverse Compton scattering. j^e prime candidates for TeV y-ray emission from this 
The well studied objects Mkn421 (redshift z = 0.030) and class of objects 
MknSOl (z = 0.034) belong to this type of TeV y-ray emitters. 

Recently, the BL Lac type objects lES 1959-1-650 (z = 0.047) The elliptical galaxy M87 (right ascension ajaoooo : 
( |Nishiyama et al. 1 999 Ah aronian et al. 2003b> and the much 12•''30'M9.4^ declination feooo.o : +12°23^28 ", redshift z ^ 
more distant H 1426-1-428 (z = 0.129) JHoran et al. 20021 0.00436) has an optical extension of 8.3^ x 6.6' dMa et al. 1998> 
with a large radio halo of 16' x 12' JCameron 1971> . 

5e«rf c#innf re^t(e.s?i' to: N. Getting, M. Tluczykont M87 contains a supermassive black hol e with a mass 

e-mail: Niels.Goetting@desy.de, Martin.Tluczykont@desy.de Mbh ~ 2-3 X 10^ GHarmS et al. 1994V The pOWer of 
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Date 


Year 


Obs. Time 








t^iJ 


r°i 
t J 


l~)ppprnhpr 


1998 


0.7 


23.6 


Tfiti 17 Tfin 

J all. LI — J dii. 


1999 


10.7 


18.4 


February 1 2 


1999 


0.7 


17.0 


March 16 -March 24 


1999 


21.7 


21.4 


Aprils- April 21 


1999 


29.4 


23.2 


May 8 -May 18 


1999 


19.9 


20.9 


June 3 


1999 


0.3 


40.5 


Total 


83.4 


21.6 



Table 1. Dates of individual HEGRA observation periods of 
M 87. Listed are observation times and mean zenith angles {ff). 
Typically, each night comprises approx. 1 - 2 h of observation 
time. 



the non-thermal jet is estimated to be as high as a few 
10'*'* ergs-' (Owen etal. 2000a I. The angle of the M87 jet 
axis to the line of sight was determined to be 30° - 35° 
( |Bicknell & Begeknan 1996 1). M 87 is located in the central re- 
gion of the Virgo cluster of galaxies, which itself is another in- 
teresting site for particle acceleration (e. g. |Volk et al. (1996)t . 

The VERITAS collaboration has targeted M 87 with the 
Whipple 10 m Cherenkov telescope in the years 2000 and 2001 
for a total time of 14 h. Positive excesses have been observed 
at low significances of 1.6 cr (2000) and 0.9 cr (2001) lead- 
ing to a 3 o- upper limit of Ny(E > 250 GeV) < 2.2 x 10"*' 
phot. cm^^s ' (ILebohec etal. 2001> . 

The HEGRA collaboration has extensively observed 
M87 in 1998 and 1999 with the stereoscopic system of 
5 imaging atmospheric Cherenkov telescopes (lACT system, 
IDaum et al. 1997> . About half of the total observation time 
(44.1 h out of 83.4 h) has been used in an earlier analysis 
( |Gotting et al. 2001) . No evidence for TeV emission was found 
in this dataset and a 3 o" upper limit on the TeV y-ray flux 
from M 87 was determined to be Ny(E > 720 GeV) < 1 .45 x 
10"'^ phot, cm^^s"'. 

In this Letter the results of the whole data set of the ex- 
tensive HEGRA M87 observations during the years 1998 and 
1999 are reported, now also applying a more sensitive analy- 
sis method. Astrophysical conclusions concerning the nature 
of the observed excess are briefly discussed. 

2. Observations and results of analysis 

M87 was observed in the years 1998 and 1999 with the 
HEGRA lACT system for a total of 83.4 h. There were no 
further observations of M 87 with the HEGRA telescopes in 
the subsequent years. The major part of the M87 data was 
taken with a 4-telescope setup. Table ^ specifies the observa- 
tion times and mean zenith angles of the individual HEGRA 
observation periods. The mean zenith angle of 21.6° can be 
converted into a mean energy threshold (defined as the peak 
detection rate for y-showers) of 730 GeV for a Crab-like spec- 
trum l |Konopelko et al. 1999> . 

Only data of good quality were considered for the analysis, 
the most critical condition being the lACT system's cosmic ray 



M 87 event selection 



stereo algorithm 


#3 




number of images per event 


>2 




shape cut on mscw 


< 1.1 




angular distance cut 


< 0.016 deg2 






background models: 




ring segment 


template 


A^ON 


716 


716 


A'oFF 


6950 


1850 


a = floN/fioFF 


0.08757 


0.325 


A^y-candidates 


107.4 ± 26.8 


114.8 ± 30.2 


significance (cr) 


4.1 


3.9 



Table 2. Cuts, event numbers, and significances for the 
HEGRA observations of M 87 resulting from the signal search 
using the ring segment and template background model, re- 
spectively (see text). 



background trigger rate not deviating more than 30 % from the 
rate expected for the current zenith angle. A total of about 5 % 
of the data was rejected due to this selection. 

All observations of M 87 were carried out in the so-called 
wobble mode targeting the object's position ("ON") as given 
in Section[2shifted by +0.5° in declination with respect to the 
center of the field of view. This observation mode allows for si- 
multaneous estimation of the background ("OFF") rate induced 
by charged cosmic rays ( Aharonian et al. 1997j . The analysis 
uses an extended OFF-region reducing the statistical error on 
the number of background events. A ring segment is chosen 
with 260° opening angle at the same radial distance to the cen- 
ter of the field of view as for the ON-source position (see also 
"Ahar onian et al . (2002byi. The width of the ring is set to the 
diameter of the ON-source area in order to provide the same 
angular acceptance for ON- resp. OFF-source events. This ring 
segment background model is similar to the usage of a set of 
control regions ( Aharonian et al. 200 It . but provides a smaller 
ratio of ON- to OFF-source solid angle areas {a - QoN/fiopp) 
and thus reduces the statistical error on the number of estimated 
background events. For a consistency check (and for a search 
for y-ray sources in the field of view, see below) the so-called 
template background model has also been used ( |RowelI (2003)| 
see also Aharonian et al. (2002c) ) . 

For the image analysis, the mirror reflectivities and pho- 
tocathode efficiencies - which degrade slowly with time - 
along with the factors converting from digitized photomul- 
tiplier signals to photoelectrons have been determined on a 
monthly basis. The shower reconstruction and the event selec- 
tion cuts have already been described in previous publications 
(e. g. [Aharonian et al. (1999)t. The stereo air show er direction 
reconstruction algorithm #3 ( Tlofmann e t al. 1999 ) and a "tight 
shape cut" (parameter mscw < 1.1) ( jKonopeUco et al. 1999, ) 
for an effective y-hadron separation have been appUed lead- 
ing to a sensitivity gain as compared to the earlier analysis of 
the HEGRA M 87 observations. The optimum angular cut was 
derived using y-ray events from the Crab nebula on the basis 
of a nearly contemporaneous data set at similar zenith angles. 
Table|2lsummarizes the event selection cuts, the resulting event 



Aharonian et al.: Is the Giant Radio Galaxy M 87 a TeV Gamma-Ray Emitter? 



3 



in 

g 380 
> 

0) 



« 360 
E 

3 
C 



340 



320 



300 



280 



260 



optimum 
angular cut 



HEGRA 

M87- 1998/1999 

• ON 

OFF 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 
squared angular distance 6^ [deg^] 

Fig. 1. Number of events vs. squared angular distance 0^ to 
the position of M87 as observed in the years 1998 and 1999 
with the HEGRA lACT system. The dots show the ON-source 
events, while the dashed histogram gives the background es- 
timate determined from a ring segment as explained in the 
text. The statistical error for the background estimate is much 
smaller than the error of the ON-source distribution. Indicated 
by the vertical dotted line is the optimum angular cut as deter- 
mined from nearly contemporaneous Crab observations at sim- 
ilar zenith angles. The significance of the M 87 excess amounts 
to 4. 1 a-. 



numbers and significances for the data set as derived using the 
ring segment and template background models. 

Figure shows the event distribution both for the ON- 
source and the OFF-source regions as a function of the squared 
angular distance of the reconstructed shower direction to the 
source position after applying all event selection cuts. The sta- 
tistical significance of the observed excess from the direction of 
M 87 (at the reference coordinates given in SectionQ) is 4. 1 cr, 
calculated using formula (17) from Li & Ma (1983)| On the ba- 
sis of the limited event statistics the excess is compatible with 
a point-like source for the HEGRA lACT system at the posi- 
tion of M 87, although extended emission cannot be excluded. 
After background subtraction the mscw values show a Gaussian 
distribution around the value of 1 .0 as expected for a y-ray pop- 
ulation (see Konopeiko et al. 1999). This test supports the hy- 
pothesis that the measured excess is a result of M 87 being a 
true TeV y-ray source. Applying different statistical tests in or- 
der to search for burstlike behaviour of M 87 no evidence for 
flux variations in the TeV energy range has been found. 

The event distribution in the field of view was used to de- 
termine the center of gravity position (CoG) of the TeV y-ray 
excess at o-jjoooo = 12'"30"54.4^ + 6.9^,,, + 1.7;,,, (5j2ooo.o = 
12°24'17"+1.7;|,^,+0.4Jyj, as shown in Figure|2] The accuracy of 
the CoG determination is limited by a systematic pointing error 
of about 25 " tPiihlhofer et al. 1997< . Within the large errors. 



HEGRA - M87 
TeV CoG position 
{la errors) 




15 00 
RIGHT ASCENSION (B1950) 



Fig. 2. Radio image of M 87 at 90 cm showing the structure of 
the M 87 halo. The center of gravity position of the TeV y-ray 
excess from the HEGRA M 87 observations is marked by the 
cross indicating the statistical 1 cr errors (radio image adapted 
from (Owen et al. (2000b)> . 



the CoG is consistent with the M 87 position, although a small 
shift of the source position cannot be ruled out. Therefore, it is 
not possible to localize a candidate TeV y-ray production site 
to particlular inner radio structures of M 87. 

In order to search for TeV y-ray sources in the relatively 
large field of view of the HEGRA lACT system (and for a con- 
sistency check with the ring segment background model), the 
template background model was used for a further analysis of 
the M 87 data. A 2° x 2° skymap of excess events determined 
for this sky region is shown in Figure |3l The excess from the 
direction of M 87 is clearly visible in the representation using 
overlapping circular bins, showing that the only significant ex- 
cess in the field of view is related to M87. The significance 
using the template model is 3.9 cr as given in Tabled The two 
background models applied in this analysis use widely differ- 
ent approaches, supporting the assumption that the M 87 excess 
does not stem from a background fluctuation. 

The observed excess can be converted into an integral flux 
of (3.3 + 0.8)% of the Crab nebula flux (only the statisti- 
cal error is given because of the low statistics). A conver- 
sion into absolute flux units using the well measured pho- 
ton flux and spectrum of the Crab nebula around 1 TeV 
(e. g. [Aharonian et al. (2000)t results in a y-ray flux of 

Ny(E > 730 GeV) = (0.96 + 0.23) x lO^'^phot. cm"^ s"^ (1) 

A spectral analysis of the data of the M87 data has 
been performed using the analysis technique described in 
[Aharonian et al. (2003a)| The data can be well described with 
a power law dN/dE ~ E'" with o- = 2.9 + 0.8,,,,, + 0.08,,,,. The 
large statistical error results from the low event statistics. 
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Fig. 3. Skymap (2° x 2° view at a 0.02° x 0.02° binning) of the 
number of excess events in the HEGRA field of view around 
the position of M 87 generated using the template background 
model. At each bin, the excess is determined from events within 
the optimum radius of 0.127° (based on Crab observations). 
Negative excesses below -15 events are shown in white color 
for clarity reasons. The solid angle resulting from the optimum 
radius is indicated by the white circle with M87 being located 
at the intersection of the white lines. 

3. Summary and Conclusions 

The radio galaxy M 87 has been observed with the HEGRA 
lACT system for a total of 83.4 h. For the first time a significant 
excess of 4. 1 cr has been detected at energies above a mean 
energy threshold of 730 GeV from a member of this class of 
objects using the imaging atmospheric Cherenkov technique. 
Note that nearly 30 years ago the radio galaxy Centaurus A 
was reported to be a TeV source observed with a non-imaging 
instrument (Grindlay et al. 1975 1. 

Due to the limited number of excess events detected by the 
HEGRA telescopes, an analysis of the spectral shape results in 
a very large statistical error, thus making it difficult to draw a 
conclusion about the origin of the TeV y-radiation. Assuming 
a power-law shape with a photon spectral index of 2.9 the inte- 
gral photon flux at the level of (3.3 + 0.8) % of the Crab nebula 
flux converts into an energy flux of 

Fy{E > 730 GeV) = (2.4 + 0.6) x 10"'^ erg cm"^ s"' . (2) 

Given the distance to M 87 of about 16 Mpc, this corresponds to 
a y-ray luminosity above 730 GeV of about lO'*' erg s"' under 
the assumption of isotropic emission. 

Several different possibilites for the origin of GeV/TeV y- 
radiation are conceivable. M 87 with its pc scale jet has re- 
cently been modeled within the Synchrotron Self Compton sce- 
nario as a BL Lac object seen at a large angle to its jet axis 
JBai & Lee 2001> . It has also been modeled using the so-called 
Synchrotron Proton Blazar model f Protheroe et al. 2002> . In 
both models, the observed flux can be accommodated. 



The large scale (kpc) jets with several knots detected at 
radio to X-ray frequencies and believed to have synchrotron 
origin due to electrons with energies up to 100 TeV is also a 
possible y-ray production site in M 87. Consequently, inverse 
Compton y-rays in the 1 - 10 TeV energy range can be ex- 
pected within reasonable model parameters. 

Moreover, y-rays could be produced in the interstellar 
medium of M 87, i. e. at larger distance scales from the center 
of this active galaxy. Both inverse Compton and hadronic inter- 
actions could generate a TeV y-ray luminosity in the observed 
range of approx. lO'*' erg s ' . If this interpretation is valid, one 
would expect a slightly extended source (2 - 3')- 

It should be noted that M 87 is also considered as a possible 
source of TeV y-rays from the hypothetical neutralino annihi- 
lation process ffialtz et al. 2000b . 

A weak signal at the centi-Crab level is at the sensitivity 
threshold for the HEGRA lACT system for observation times 
of the order of 100 h. Therefore, a deep investigation and pos- 
sible confirmation with a spectral analysis of the M 87 excess 
should be subject of the next generation Cherenkov telescope 
projects like H E S S, MAGIC and VERITAS, which provide 
increased sensitivity together with a lower energy threshold. 
Due to the proximity to M 87 (16 Mpc) and to the increased ac- 
curacy of these observations (a fraction of an arc minute), these 
measurements may allow the location of the y-ray production 
site in M 87 to be more accurately determined thus greatly ad- 
vancing our understanding of its TeV y-radiation. 
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